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INTRODUCTION 

Irrigation  of  agricultural  lands  is  one  of  the  greatest  advancements 
that  man  has  achieved  in  his  ever  increasing  battle  against  the  unfavor- 
able conditions  that  stand  in  the  way  of  good  and  enormous  agricultural 
production. 

In  irrigation,  like  many  other  fields  of  human  endeavor,  the  struggle 
to  move  the  system  toward  perfection  has  been  the  concern  of  many.  Engineers 
working  in  the  area  have  been  forced,  by  the  constraints  that  most  irriga- 
tion systems  provide,  to  redouble  their  efforts  toward  a  higher  efficiency 
in  the  use  of  the  single  most  important  resource  required,  for  the  survival 
of  the  system,  namely  water. 

Water  losses  due  to  evaporation  and  seepage  in  open-ditch  surface  ir- 
rigation systems  have,  for  a  long  time,  been  matters  of  concern.  The 
advent  of  pipe  irrigation  systems  have  proved  a  remedy  to  those  problems 
but  not  without  presenting  some  of  their  own.  In  furrow  irrigation,  for 
instance,  gated  pipes  replace  the  ditches  and  siphon  tubes  for  supplying 
water  to  furrows.  However,  the  gated  pipes  present  a  problem  of  non-uniform 
flow  through  the  various  openings  of  the  pipe.  The  ramifications  of  the 
problem  are  tantamount  to  an  inefficient  use  of  the  water.  If  we  take, 
for  example,  the  condition  in  which  more  water  comes  out  of  the  inlet  end 
gates  than  from  the  closed  end  gates,  the  tendency  would  be  either  to 
irrigate  adequately  at  the  beginning  and  under-irrigate  towards  the  end 
furrows  along  the  pipe  or  to  over-irrigate  at  the  beginning  furrows  and 
irrigate  adequately  at  the  end  furrows.  The  phenomenon  reverses  when  there 
is  more- water  coming  out  of  the  end  furrows  than  at  the  beginning.   In 


either  case  the  consequences  are  undesirable.  At  present,  farmers  manually 
adjust  the  gates  in  an  attempt  to  obtain  uniform  flow;  an  exercise  that 
requires  skill  and  good  judgement  and  also  one  that  shall  always  produce 
a  result  which  is,  at  best,  only  approximate.  Moreover,  in  places  where 
labor-saving  is  an  element  of  major  concern  to  the  economic  well  being  of 
the  system,  there  is  a  desire  to  automate  the  system.  This,  at  present, 
is  practical  except  for  the  problem  of  attaining  uniform  flow  out  of  the 
gated  pipes.  Since  with  an  automated  system  the  aim  is  to  minimize  labor 
input,  it  shall  be  necessary  to  find  a  means  of  obtaining  a  uniform  flow 
through  the  gates  without  the  need  for  making  adjustments. 

The  solution  to  this  problem  requires  a  careful  examination  of  the 
factors  that  influence  flow  out  of  the  various  orifices  of  the  pipe. 
The  two  factors  that  are  known  to  have  a  direct  effect  upon  the  flow  are 
the  static  head  at  the  orifice  and  the  discharge  coefficient  of  the  ori- 
fice. The  possibility  exists  that  this  latter  factor,  the  discharge  co- 
efficient, can  in  turn  be  influenced  by  several  other  factors.  In  fact, 
investigations  and  analysis  of  the  problem  made  to  date  have  indicated 
the  presence  of  such  factors.  Two  of  them  that  have  received  an  almost 
unanimous  recognition  are  the  velocity  of  water  inside  the  pipe  as  it 
approaches  the  orifice  and  the  static  head  at  the  orifice.  However,  the 
nature  of  the  relationship  between  these  factors  and  the  coefficient  has 
not  been  clearly  understood  or  established. 

The  purpose  of  this  study  was  to  carry  out  an  investigation  to 
determine  the  relationship  between  the  coefficient  of  discharge  for 
orifices  in  a  gated  pipeline  and  the  two  factors,  velocity  of  water  ap- 
proaching the  orifice  and  static  head  at  the  orifice. 


REVIEW  OF  LITERATURE 

The  problem  of  obtaining  a  uniform  flow  through  several  openings  on 
the  side  of  a  pipe  is  not  unique  to  gated  irrigation  pipes.  It  has  been 
of  concern  to  engineers  dealing  with  manifolds  in  different  systems  well 
before  the  development  of  gated  irrigation  pipes.  These,  among  others, 
include  the  pipe  burner  for  gaseous  fuels,  manifold  of  the  radiant  fire 
headers  for  air  heaters  and  distribution  pipes  in  water  filtering  systems, 

Some  of  the  work  was  more  or  less  of  a  trial  and  error  nature  in 
which  the  goal  was  to  obtain  uniform  discharge  through  orifices  of  the 
manifold.  Attention  in  earlier  works  had  not  been  paid  specifically  to 
the  factors  on  which  the  uniformity  of  discharge  depends. 

As  early  as  the  middle  of  the  nineteenth  century,  Francis  (1865) 
gave  rules  for  the  relation  of  orifice  to  pipe  size  in  order  to  accom- 
plish uniform  discharge  in  the  specific  case  of  perforated  pipes  in  a 
fire  protection  system.  However,  he  did  not  analyze  the  various  factors 
involved  in  the  variation  of  the  discharge  from  one  orifice  to  another. 
Later,  Jenks  (1922)  published  more  elaborate  rules  of  a  similar  kind. 
These  seemed  to  be  confirmed  by  the  results  of  studies  on  pressure  dis- 
tribution conducted  by  Elms  (1927). 

Enger  and  Levy  (1929)  published  results  of  experimental  studies 
which  indicated  that  the  Bernoulli  equation  could  be  used  to  establish 
the  pressure  distribution  at  various  points  in  a  pipe.  The  results 
further  revealed  the  effect  of  water  velocity  inside  the  pipe  on  the 
coefficient  of  discharge  through  orifices.  In  their  analysis,  they 
negelected  the  pipe  friction  and  hence,  determined  that  pressure  at  any 


point  on  a  slot  is  equal  to  the  head  at  the  end  of  the  slot  minus  the 
velocity  head  at  the  given  point.  This  they  concluded  from  the  follow- 
ing equation: 

v  =Y2g  (H-y)  (1) 

where 

v  =  velocity  in  the  pipe  at  any  section  in  question 

H  =  head  at  the  end  of  the  slot 

y  =  pressure  at  any  point  in  question 

g  =  acceleration  due  to  gravity 

They  derived  the  above  equation  from  the  assumption  of  zero  friction,  the 

impulse  momentum  principles  and  the  boundary  condition  that  at  the  end  of 

the  slot,  velocity  inside  the  pipe  was  zero. 

Assuming  a  constant  discharge  coefficient,  Enger  and  Levy  derived 

another  expression,  this  time  for  pressure: 

y   =  H  vers  (lf  _  2Ce  (L_X)  )  (2) 

where 

C  =  coefficient  of  discharge  of  the  slot 

e  =  slot  width 

L  =  slot  length 

X  =  a  distance  from  the  beginning  of  the  slot  to  a  point  in  question 
However,  they  later  conceded  to  the  fact  that  the  discharge  coefficient 
was  not  constant  but  rather  decreased  with  an  increase  of  the  velocity 
in  the  pipe  in  accordance  with  the  empirical  formula: 

v2 

y  ~   27T  (3) 

C  =( £fl-)  Ce 

y 


where 

C  =  coefficient  of  discharge  on  any  opening 

Ce  =  coefficient  of  discharge  of  the  last  opening 
Computed  values  for  the  discharge  coefficient  using  Equation  3  were 
found  to  be  in  close  agreement  with  experimental  values  obtained  in 
other  studies.  It  is  worth  noting  here  that  Enger  and  Levy  did  address 
themselves  to  the  relation  of  orifice  coefficient  with  both  the  static 
and  velocity  heads  inside  the  pipe.  Nevertheless,  it  should  be  pointed 
out  that  in  their  entire  analysis  they  did  not  take  into  account  the 
rate  of  recovery  of  pressure  due  to  deceleration  of  flow  in  the  axial 
direction  which  is  due  to  the  redirection  in  fluid  momentum  as  a  result 
of  the  discharge  through  the  various  orifices.  Further,  they  didn't 
consider  the  effect  of  friction  on  the  discharge  rate. 

Wills  (1931)  had  performed  some  investigations  on  gas  flow  through 
burner  manifolds.  He  also  neglected  friction  losses  and  regarded  the 
manifold  as  an  infinite  reservoir.  Thus,  the  cross-sectional  area  of  the 
manifold  is  a  lot  larger  than  the  cross-sectional  areas  of  all  the  ori- 
fices combined.  Hence  from  these  assumptions,  uniformity  of  flow  was 
eminent. 

Keller  (1949), in  a  comprehensive  paper,  dealt  with  the  case  of  a 
manifold  supplying  fluid  to  a  set  of  parallel  pipes  or  ducts,  or  dis- 
charging through  numerous  orifices  distributed  along  a  pipe.  He  did  not 
assume  any  particular  state  of  flow,  neither  laminar  nor  turbulent.  He 
took  the  friction  factor  to  be  constant  along  the  entire  length  of  the 
manifold  irrespective  of  the  flow  rate.  He  also  regarded  as  perfect  the 
pressure  recovery  due  to  the  axial  deceleration  of  the  flow.  In  addition, 


he  neglected  pressure  losses  due  to  branching.     In  his  calculations,  Keller 
determined  the  variation  in  width  of  a  longitudinal   slot  in  a  pipe  which 
would  provide  uniform  outflow.     He  also  determined  the  required  change  in 
cross-sectional   area  of  a  pipe  essential   for  the  provision  of  uniform  out- 
flow through  slots  of  constant  width.     This  latter  aspect  was  treated  by 
W.  M.   Dow  as  will   be  mentioned  later. 

The  equation  obtained  by  Keller,  for  the  variation  of  discharge  along 
a  pipe  length  of  constant  cross-sectional   area, is: 

vi  ■  swim  '  C05  (KR  S/L)  (4) 

where 

S  =  distance  from  the  dead  end  of  the  manifold 

VQ  =  inlet  velocity 

Vj  =  discharge  velocity  at  distance  S  from  the  dead  end  of  the  manifold 

K  =  coefficient  of  discharge  of  the  holes  or  slot 

L  =  length  of  the  manifold 

R  =  area  ratio  =  sum  °f  areas  of  all  discharge  openings 

cross-sectional   area  of  manifold 

It  can  be  noted  that  the  above  equation  is  a  form  of  the  combination  of 

Equations  1  and  2  given  earlier  as  determined  by  Enger  and  Levy. 

Dow  (1950),  in  his  theoretical  analysis  of  flow  through  a  perforated 

pipe  with  a  closed  end,  assumed  that  it  was  possible  to  obtain  a  uniform 

discharge  throughout  the  length  of  the  pipe  when  the  spacing  and  diameter 

of  the  orifices  in  the  pipe  wall  were  uniform.  He  also  assumed  a  perfect 

rate  of  pressure  recovery.  His  analysis,  however,  was  restricted  to 

finding  the  cross-sectional  area  of  the  pipe  as  a  function  of  the  distance 

along  its  axis.  An  implicit  assumption  that  can  be  gathered  from  Dow's 


analysis  is  that  if  uniform  discharge  is  to  be  attained  with  constant 
orifice  area  and  spacing,  then  of  necessity,  the  cross-sectional  area  of 
the  manifold  has  to  vary.  Hence,  he  determined  the  area  as  a  function 
of  the  axial  length.  The  results  of  his  analysis  were  verified  by  the 
photographs  of  flame  heights  from  gas  burners  both  before  and  after  modi- 
fication. This  was,  however,  restricted  to  relatively  short  pipes. 

The  manifold  problem  was  considered  in  a  relatively  more  detailed 
form  by  Van  Der  Hegge  Zijnen  (1951).  One  of  the  typical  cases  of  the 
problem  he  discussed  was  that  of  a  turbulent  flow  inside  the  manifold 
and  turbulent  discharge  through  the  orifices.  The  assumptions  he  made 
were  those  of  isothermal  conditions,  incompressibility  of  the  fluid,  con- 
tinuity of  discharge  over  the  manifold  length  and  a  continuous  cross- 
sectional  area  of  the  manifold.  He  further  assumed  a  constant  coefficient 
of  pressure  recovery  and  constant  coefficient  of  discharge. 

Van  Der  Hegge  Zijnen,  in  his  analysis,  started  with  the  requirements 
that  to  obtain  uniform  discharge  per  unit  length  of  the  manifold,  the 
orifice  spacings,  /n,  could  be  determined  if  the  diameters  of  both  the 
pipe  and  the  orifices  were  kept  uniform,  or  the  orifice  diameter  could 
be  determined  with  the  pipe  diameter  and  orifice  spacings  made  uniform 
or  else  the  pipe  diameter  at  any  point  could  be  determined  while  the  ori- 
fice diameters  and  their  spacings  were  kept  constant.  He  assumed  that 
the  initial  discharge  into  the  pipe  and  pressure  head  at  inlet  were  always 
known  or  could  be  prescribed.  Among  the  conditions  he  considered  were 
pressure  recovery  due  to  deceleration  of  flow  as  water  exits  through  the 
orifices  which  he  assumed  counteracts  the  fluid  friction  and  the  bending 
of  flow  from  the  axial  direction  of  the  manifold  to  the  lateral  direction 


8 


of  the  discharge.  The  latter  condition,  he  believed,  would  make  fluid 
resistance  accompany  the  discharge  through  the  orifices.  He  further  sug- 
gested that  the  magnitude  of  the  fluid  resistance  was  determined  by  the 
velocity  of  flow  just  upstream  from  the  orifice  in  question.  He  started 
his  analysis  with  the  equation  of  motion  for  the  flow  in  a  manifold,  which 
was  given  as: 

A     (p+K-%PV2)  =-A-%PV2/D  (5) 

dx 

where 

p  =  pressure  just  upstream  from  an  orifice 

K  =  rate  of  pressure  recovery 

P  =  water  density 

V  =  velocity  upstream  from  the  same  orifice 

A  =  coefficient  of  fluid  friction 

D  =  diameter  of  the  pipe 

For  laminar  flow  in  a  pipe  discharging  through  orifices,  Van  Der 

Hegge  Zijnen  set  out  to  determine  either  the  spacing  or  the  diameter  of 

the  orifice  per  unit  length,  given  that  one  of  them,  was  uniform  and  the 

pipe  diameter  was  the  same  everywhere.      In  order  to  obtain  uniform  flow, 

the  following  condition  had  to  be  satisfied: 

q  =     J    d2  a[PQ  -  %PV2   (fi-     -  2K  ft)     J       (1   -  £)lh{Z/P)h  (6) 

o 

where 

q  =  discharge  through  each  of  the  orifices 

P  =  the  pressure  at  the  inlet 

Re    =  the  Reynold's  number  at  the  inlet 
o 

a  =  a  proportionality  constant 
x  =  distance  from  the  inlet  end 
L  =  the  total  length  of  the  pipe 


The  discharge  q  through  the  various  orifices  must  sum  up  to  the  discharge 

Q0  into  the  pipe  at  the  inlet.     Hence     qQ  =  Q0/nL  in  which  n  is  the 
+ 

number  of  orifices  per  unit  length.     Equation  6  can  then  be  used  to  deter- 
mine the  spacing,  -  ,  when  both  the  pipe  and  the  orifice  diameters  are 
uniform.     On  the  other  hand,   the  same  equation  can  be  used  to  determine 
the  various  sizes  of  the  orifices  if  the  spacing  and  pipe  diameter  are 
made  uniform.     In  the  case  of  turbulent  flow,   the  following  equation  had 
to  be  satisfied  for  uniform  flow: 

Po  +¥n,o2    [2  JSf     (1  -  £)  +  0.158  ReQ     ~\  £      [ft-jL)  f  -  1]  ]  = 

Q2  f     [(  7r     d2nLa)2  2/p-j  (7) 

The  other  alternative  that  Van  Der  Hegge  Zijnen  considered  was  that 

of  maintaining  a  uniform  spacing  and  equal   orifice  sizes  while  the  pipe 

diameter  changes  per  unit  length.     He  developed  the  following  equation: 

Re  "^ 
A_    =  i  +  0.158     ~-^~    [l-(l-x/L)3/8]  8/3       (1-x/L)  (8) 

where 

A  =  cross-sectional  area  at  any  point 

AQ  =  cross-sectional  area  at  the  inlet 
Gladding  (1940)  discussed  the  loss  of  head  in  a  pipe  of  uniform 
diameter.  Using  a  simple  analysis  made  possible  by  the  assumption  that 
it  was  possible  to  obtain  an  equal  discharge  through  evenly  spaced  out- 
lets on  the  side  of  a  pipe,  he  showed  the  loss  of  head  can  be  given  by 
the  expression: 


h  =  FKLQm  (9) 


where 

h  =  loss  of  head 
K  =  a  constant 
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m  =  a  constant 

L  =  length  of  the  pipe 

Q  =  total  discharge  into  the  pipe 

F  =  a  quantity  which  itself  is  given  by 

F  =  ZN  '/Nm    in  which  N  is  the  number  of  outlets 
Howl  and  (1953)  presented  a  paper  to  the  3rd  Midwestern  Conference 
on  Fluid  Mechanics  in  which  he  discussed  the  design  of  perforated  pipe 
for  uniformity  of  discharge.  He  established  from  Bernoulli's  equation 
the  variation  of  pressure  head  along  a  manifold  pipe  with  closed  end 

when  uniformity  of  discharge  was  attained.  He  gave  the  expression  as: 

Vl       2 
h  =  ho  "  W     (X/V  +  hL  (10) 

where 

h  =  pressure  head  at  any  point  along  the  pipe 

hQ  =  pressure  head  at  the  closed  end  of  the  pipe 

V  =  velocity  at  entrance  to  the  pipe 

h.  =  head  loss  due  to  pipe  friction  in  the  portion  of  the  pipe  between 
the  point  under  consideration  and  the  closed  end 

X  =  distance  from  the  point  to  the  closed  end 

X,=  length  of  the  pipe 
Howland  further  showed  that  the  friction  loss,  h,  ,  could  be  expressed  as: 

hL  =  0^   S  dx  =  o/X  Sdx  ("> 

where 


s  =dh 


dx  "  Sl  ^V" 


Sj  =  head  loss  per  unit  length  at  the  entrance 

S  =  head  loss  per  unit  length  at  point  under  consideration 

V  =  velocity  at  point  under  construction 

n  =  a  constant  varying  between  1.75  and  2.00 
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The  final  form  of  the  total  head  loss  equation  was: 

V,2  f,X, 


h"ho    :W      t  (Mtd    (x/xi)3-  (x/xi)2]  (12) 


where 


fi   =  friction  factor 
D     =  pipe  diameter 

Howland  gave  consideration  to  the  change  of  orifice  coefficients 
with  a  change  in  velocity.     He  developed  an  empirical   curve  that  estab- 
lished the  orifice  coefficient  as  a  function  of  pipe  velocity. 

Hansen  (1954)  wrote  a  paper  on  the  determination  of  water  flow  from 
gated  pipe  in  which  he  indicated  that  uniformity  of  flow  was  obtained 
when  he  used  two  20-foot,  4-inch  gated  pipes  connected  at  the  middle  by 
a  5-foot  tee  section  and  placed  on  a  1  in  300  slope.     He  used  22-inch  and 
36-inch  gate  spacings  to  obtain  data  to  establish  discharge-opening-head 
relationships. 

Tovey  (1959)   followed  the  recommendation  of  Hansen  using  similar 
4-inch  gated  pipe  and  attempted  to  duplicate  his  results  but  found  that 
uniform  flow  could  not  be  obtained  for  0.55-inch  gate  openings     where 
the  average  head  loss  per  gate  exceeded  0.1  foot.     In  a  similar  manner 
he  also  found  that  with  0.95-inch  and  1.7-inch  gate  openings,  uniform 
flow  could  not  be  obtained  when  the  average  head  loss  per  gate  exceeded 
0.15  foot  and  0.17  foot,   respectively. 

Spomer  (1961)   also  used  the  Hansen  slope  recommendation  for  a  6-inch 
pipe.     He  found  that  instead  of  uniform  flow,  a  rising  slope  of  1  in  300 
gave  a  decrease  in  discharge  from  the  gates  as  one  approaches  the  dead 
end.     The  reverse  occurred  when  the  pipe  was  on  a  falling  slope. 
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Ramirez-Guzman  and  Manges   (1971),   in  a  paper  presented  to  the  annual 
meeting  of  the  American  Society  of  Agricultural   Engineers,  discussed  uni- 
form flow  from  orifices  in  irrigation  pipe.     They  used  a  method  of  analysis 
similar  to  that  of  Howland  mentioned  earlier.     They  also  obtained  an 
equation  similar  to  Equation  12  derived  by  Howland.     The  major  difference 
between  the  two  analyses  was  the  difference  in  the  initial   conditions  of 
the  pipe.     Howland  considered  the  pipe  to  be  horizontal  while  Ramirez-Guzman 
and  Manges  took  the  pipe  to  be  sloping.     Hence,  the  elevation  term  appeared 
in  the  latter 's  derivation.     Further,  while  Ramirez-Guzman  and  Manges  used 
the  Hazen-Williams  equation  to  express  the  energy  slope  as  a  function  of 
the  velocity,   the  friction  coefficient  and  the  hydraulic  radius  of  the 
pipe,  Howland  expressed  the  same  quantity  (i.e.   the  slope)   as  a  function 
of  pipe  friction,  the  pipe  diameter  and  the  flow  velocity  inside  the  pipe. 
The  following  two  equations  express  slope  as   used  by  Ramirez-Guzman  and 

Manges,  and  Howland,  respectively: 

Vj  1.85         . 

hi    =  L  i     oi  or        -J 


'1  ~L1.318CH     J  R1.17  .  (13) 

where 

h,  =  energy  slope  at  pipe  entrance 

Vj  =  velocity  at  the  entrance  to  the  pipe 

C^  =  Hazen  Williams  friction  coefficient 

R  =  hydraulic  radius 

2 

(14) 


1    2g  D 

where 

fi 

=  pipe  friction 

vi 

=  velocity  at  entrance 

D 

=  pipe  diameter 
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The  final  form  of  the  equation  derived  by  Ramirez-Guzman  and  Manges 

was:       o 

V/       o       V     1.85    X  2.85   qy 

h  "  ho  "  2T  ^   +  ^318C^    ^OT  W*l>     +  f  Xl   (15) 


where 

h  =  pressure  head  at  any  point  in  the  pipe 

h  =  pressure  head  at  dead  end 

X,  =  distance  of  the  inlet  end  from  the  dead  end 

X  =  distance  from  the  dead  end  to  the  point  in  question 
Using  the  above  equation,  Ramirez-Guzman  and  Manges  calculated  the 
head  at  various  points  along  the  pipe  and  thence,  determined  the  discharge 
out  of  each  orifice  assuming  a  unity  coefficient  of  discharge.  They  then 
went  back  to  correct  for  the  assumption  of  unity  coefficient  of  discharge. 
This  time  around,  they  took  the  coefficient  to  be  the  ratio  of  the  combined 
discharges  out  of  the  orifices  to  the  total  discharge  that  entered  the 
pipe.  Based  on  this,  they  adjusted  their  calculations  for  the  discharge 
through  each  of  the  orifices.  Upon  comparison  with  experimentally  deter- 
mined data,  as  observed  by  Spomer,  some  agreements  were  achieved  with 
variation  of  up  to  6.5%.  They,  therefore,  concluded  that  the  equation 
they  derived  was  accurate  enough  for  field  use. 

Even  though  in  their  calculations  Ramirez-Guzman  and  Manges  had 
assumed  that  the  orifice  coefficients  were  the  same,  they  nonetheless 
conceded  that,  due  to  the  consistent  overestimation  of  the  orifice  dis- 
charge at  the  inlet  end  and  an  underestimation  near  the  dead  end,  there 
was  an  inverse  variation  in  the  magnitude  of  the  discharge  coefficient 
with  distance  from  the  dead  end. 
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Chu  and  Moe  (1971)  followed  up  the  analysis  made  by  Ramirez-Guzman 
and  Manges.  They  introduced  a  new  element  by  taking  into  consideration 
the  variation  in  the  coefficient  of  discharge  coefficient  with  velocity. 
It  should  be  mentioned,  however,  that  their  initial  approach  differs  from 
that  of  Ramirez-Guzman  and  Manges  even  though  it  was  obviously  based  on 
the  Tatter's  analysis. 

Chu  and  Moe  started  with  Bernoulli's  energy  equation  expressed  be- 
tween two  points  at  the  centers  of  two  adjacent  orifices.  Thus  they  wrote: 

hm  *  W1    +  zm  ■  *,  ♦  tf  ♦  Vl      (16) 

whe  re 

Z  =  elevation  at  a  point 

i  =  1.  2 ,  N  =  the  subscript  identifying  the  outlet 

N  =  the  total  number  of  outlets 

h^  =  head  loss  due  to  friction  and  branching  of  flow 

Flow  through  individual  orifices  can  be  expressed  as: 

. (17) 

q  =  ca  i   2gh 

where 

q  =  discharge  out  of  an  orifice 

c  =  orifice  coefficient  of  discharge 

a  =  area  of  the  orifice 

h  =  head  at  each  orifice 

g  =  acceleration  due  to  gravity 
Chu  and  Moe  pointed  out  that,  in  order  to  have  uniform  flow  from  the  various 
orifices  at  constant  spacings  and  equal  sizes,  the  head  at  each  outlet 
should  be  the  same,  if  the  orifice  coefficient  remains  constant.  Thus: 
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hi+1  =  h1  (18) 


where 

-;     =  0  to  N 
The  energy  equation  is  thus  reduced  to: 


vni  vi2  <19) 

+  Z_. .:       =     J-      +  Z.   +  h. 


2g  i+1  2g  i         L 

They  further  suggested: 

2  2 

hL  ■  hf  +  f  t^p-  +  zH  (20) 

where 

hf  =  the  frictional  loss  of  head  of  pipe 
K  =  coefficient  of  energy  loss  of  head 
It  is  worth  noting  here  that  Chu  and  Moe  took  into  consideration  the 
loss  that  resulted  from  the  branching  of  flow.  In  the  last  part  of  their 
analysis,  they  removed  the  assumption  that  the  orifice  coefficients  are 
the  same  for  all  of  the  orifices. 

In  their  analysis,  Chu  and  Moe  developed  a  factor  which  is  a  function 
of  the  number  of  outlets  on  a  lateral.  They  expressed  the  factor,  F(N;  CF, 
K)  as: 

F  (N;  CF,  K)  =  ^j-  |  N3  -  (1+|)  N2  +  (l  +  ff)  N-  {   +  (F  "  1)   (2l) 

where 

CF  =  k   F,2/Sf 
1    1 

F,  =  Froude  number  with  respect  to  1st  section  of  the  lateral  counting 

from  the  dead  end 

Sf  =  the  slope  of  the  energy  grade  line  for  the  same  section 

1 
F   =  Christianson's  factor 
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Also: 


F   (N;   CF,    K)   =   (ZN  -   Z1)/SfN   (NL) 


(22) 


where 


S^N  =  head  loss  due  to  friction  in  the  corresponding  main  pipeline 


'fN 


Sf1   L(N)m+1 

JL    +    JL    +     (m-D: 

m+1  2N  ,M2 

6N 


(22a) 


m      =  constant    =1.9 

A  table  of  values  for  F  has  been  generated  by  Christianson   (1942) 
for  different  values  of  N  and  K. 

Substituting  Equation  22a  into  Equation  21  and  simplifying  yields  the 
following: 

rrri/      M'-~m  i    -r,  i  i 

(23) 


F  (N;   CF,   K)     •  S*    N  ""  -  CF  N1_m       +     1 


1 


m+1  2N 


Chu  and  Moe  finally  concluded  with  an  expression  for  the  slope  of 
a  pipe,  S    ' ,  where 

V  =  FT    [F(jI;  CF'  K)  Sf  jl     "  F(M>I+1J  CF>  K)  Sf(j-l)l+l^     ^24) 


where 


S  '  =  slope  of  the  j   pipe  counting  from  the  dead  end 
j   =  the  number  that  identifies  the  gated  pipe 

sP  =  y/(1-y2)Js 

X.  L. 

S   =  slope  of  the  j   pipe  counting  from  the  dead  end,  with  respect 
to  ground  surface 
The  slope, S',  was  also  regarded  as  the  average  slope  on  which  the  pipe 
should  be  installed. 

In  a  paper  published  four  years  later  by  the  Journal  of  Mechanical 
Engineering  Science,  Bailey  (1975)  discussed  and  analyzed  the  nature  of 
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"fluid  flow  in  perforated  pipes."  Bailey  pointed  out  two  factors  as 
being  the  major  causes  of  variation  of  flow  along  a  perforated  pipe  or 
a  manifold,  namely  friction  and  decrease  in  fluid  momentum.  While  a  lot 
of  work  had  been  accomplished  in  the  field  of  friction  losses  in  pipes, 
he  stated,  there  was  little,  if  anything  done  or  known  about  the  nature 
in  which  the  decrease  in  fluid  momentum  affected  the  fluid  pressure  down- 
stream of  an  outlet.  Elaborating  on  the  point,  he  maintained  that  the 
decrease  in  fluid  momentum  resulted  in  an  inevitable  decrease  in  the 
fluid  velocity  inside  the  pipe  after  passing  each  discharging  orifice. 
This,  he  said,  resulted  in  a  corresponding  increase  of  the  fluid  pressure, 
a  phenomenon  variously  described  as  the  "diffusion,"  "inertia"  or  the 
"static  regain"  effect. 

In  his  analysis,  Bailey  suggested  that  the  component  of  the  velocity 
of  the  emerging  fluid,  which  was  perpendicular  to  the  plane  of  the  ori- 
fice, was  derived  from  the  excess  static  pressure  in  the  pipe.  The  fluid 
was  assumed  to  be  coming  out  at  an  angle  to  a  perpendicular  line  which, 
in  turn,  is  at  right  angles  to  the  plane  of  the  orifice.  The  velocity 

component,  Vy,  was ,  therefore,  expressed  accordingly  as: 

_  2(P1  +  P2)   „ 

y     L      j 


where 


V     =  y-component  of  the  emerging  fluid  velocity 

P.   =  pressure  inside  the  pipe  in  the  pipe  section  preceding  the 

orifice  in  question 
?2  -  pressure  in  the  pipe  after  passing  the  same  orifice 
p    =  fluid  density    which  is  constant 
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It  can  be  noted  that  Bailey  had  assumed,  by  implication,  that  the  pres- 
sure, P,  at  the  orifice  is  the  average  of  the  pressure  upstream  and  the 
one  downstream  of  the  orifice  (i.e.  P  =  (P,  +  P?)/2).  Equation  25  above 
can  hence  be  reduced  to  the  one  usually  obtained  by  the  Terricelli's 
principle  (i.e.  Vy  =  2gh,  where  h  =  P/pg).  Bailey  then  described  the 

flow  out  of  each  orifice  as: 

a*  P,  +  Pa  , 

^  =  T  Cd  [  "V1  ]  (26) 

where 

q  =  the  discharge  out  of  the  orifice 

C  ,  =  discharge  coefficient 

d  =  the  orifice  diameter 
As  a  result  of  this  discharge,  he  explained  that  the  momentum  balance 
parallel  to  the  pipe  axis  can  be  established  as: 

P2  "  p!  =  P  (V22  -  V22)  (27) 

where 

Vj  =  velocity  in  the  pipe  prior  to  the  orifice  of  interest 
V^  =  velocity  in  the  pipe  after  passing  the  same  orifice 
The  relation,  however,  was  established  upon  the  assumption  that  the 
emerging  fluid  had  lost  all  of  its  axial  momentum.  That  is,  it  was  emerg- 
ing at  90  to  the  axis  of  the  pipe.  However,  since  this  was  not  so  in 
reality,  a  correction  factor  was  introduced  into  the  above  expression. 
Thus: 

P2  "  Pi  "  -J   P  (Vx2  -  V22)  (28) 

?2  ~   P}  is  the  pressure  increase  across  the  orifice  which  has  been 
referred  to  as  the  'static  regain'.  The  coefficient  C  introduced  in  the 
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above  equation  is,  therefore,  defined  as  the    'coefficient  of  static  re- 
gain'.    Rearranging  the  equation,   he  showed  the  coefficient  can  be  given 
by: 

C     =     P2  -  Pl  (29) 

r       MVj2-  v22) 

Further  to  these    analyses,  Bailey  had  carried  out  an  experiment  using 
two  lengths  of  0.2m  diameter  rigid     PVC       tubing,  each  of  length  1.2m. 
The  two  were  connected  with  flanges.     At  0.25m  downstream  from  the  flange 
a  25mm  orifice,  punched  in  a  0.064mm  thick  polyethelene  sheet  was  mounted 
beneath  a  100mm  diameter  hole  made  on  the  wall   of  the     PVC       tubing.     In 
a  plane  perpendicular  to  the  axis  of  the  orifice,  static  pressure  tap- 
pings were  made  on  the  surface  along  one  side  of  the  test  section.     Air 
was  used  as  the  system's  fluid.     A  pitot  probe  was  used  to  determine  the 
dynamic  head  of  the  emerging  air  with  the  aid  of  a  collecting  chamber. 
A  pitot  traverse  was  used  to  determine  the  dynamic  head  inside  the  pipe. 

As  a  result  of  his  experiment,  Bailey  came  up  with  an  empirical 
equation  expressing  the  orifice  discharge  coefficient  as: 

Cd  =  0.62  +  0.070  3  +  0.08832  (30  ) 

where 

P     +  P 
B  '  log     [log  [(1  +     '  ,     2)]  ] 

V 

Since  the  distributor  he  considered  in  his  paper  consisted  of  a  large 
number  of  diametrically  opposed  pairs  of  perforations,  Bailey  found  it 
impractical   to  measure  the  changes  in  velocity  after  passing  each  hole, 
because  of  the  small   amounts  of  air  coming  out  of  each.     Hence,  he  deter- 
mined the  static  pressure  profile  experimentally  and  developed  a  model 
with  the  aid  of  computer  that  gave  the  coefficient  of  static  regain  as: 
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V 
C  =  0.78  +  *log  [v    n  ]  (31) 

<J>  can  be  expressed  as: 

4>  =  0.284  +  0.098  log  |  (32) 

where 

d  =  diameter  of  an  orifice 

D  =  diameter  of  the  perforated  pipe 
In  addition,  Bailey  was  able  to  measure,  by  means  of  an  anemometer, 
the  angle  of  discharge  for  each  orifice.  He  subsequently  established  the 
following  empirical  relationship: 


V  ?p   . 
^P1  +  P2>  (33) 

where 


e  =  9  arctan  (^ — — ^ 


6  =  0.71  +  0.0043d,  the  angle  that  the  fluid  made  with  the  axis  of 

the  orifice 
d  =  diameter  of  the  orifice 
Vi=  approach  velocity  to  the  orifice 
P,  =  pressure  in  the  section  just  before  the  orifice 
?2  =  pressure  just  after  the  orifice  at  the  section  immediately  down- 
stream from  it 
In  conclusion,  Bailey  suggested  that  the  use  of  these  coefficients 
will  enable  predictions  to  be  made  of  both  the  static  pressure  and  fluid 
discharge  variations  along  uniformly  perforated  distributors.  He  main- 
tained that  such  predictions  had  agreed  with  experimental  observations. 


INVESTIGATION 
Objectives 

Analyses  presented  earlier  and  studies  by  other  researchers  have 
suggested  the  presence  of  some  effect  of  an  approach  velocity  towards 
an  orifice  and  of  pressure  head  at  the  orifice  upon  its  discharge  coef- 
ficient. The  review  of  literature  revealed  the  existence  of  only  one 
empirical  formula  relating  these  quantities.  The  formula,  given  earlier 
as  Equation  3,  cannot  be  used  without  assuming  some  value  for  the  dis- 
charge coefficient  at  the  dead  end.  However,  since  the  value  at  the  dead 
end  is  not  constant,  the  criterion  for  choosing  it  shall  have  to  be  ar- 
bitrary. Hence,  the  accuracy  of  any  prediction  of  the  discharge  coef- 
ficient at  points  other  than  the  dead  end  shall  have  to  depend  upon  the 
accuracy  of  the  assumed  value  at  the  dead  end.  Therefore  this  investi- 
gation was  undertaken  with  the  following  objectives. 

1.  To  determine  the  relationship  between  orifice  discharge  coef- 
ficient and  two  parameters;  approach  velocity  and  static  head. 

2.  To  determine  if  the  results  of  studies  by  Enger  and  Levy  on 
2-inch  diameter  pipes  can  be  applied  to  larger  pipe  sizes,  such  as 
those  used  for  gated  pipe. 

3.  To  determine  if  equations  obtained  as  a  positive  achievement 
of  objective  1  fit  other  published  data  of  flow  from  orifices  of 
a  manifold  pipe. 

4.  To  determine  a  procedure  for  designing  a  system  that  shall  pro- 
vide uniform  flow  out  of  the  various  outlets  of  a  gated  pipe  system. 
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Theory 

Whenever  there  is  an  opening  in  the  wall  of  a  conduit  or  any  con- 
tainer such  that  a  fluid  issues  through  the  entire  area  of  the  opening 
under  pressure,   the  opening  is  referred  to  as  an  orifice. 

The  conservation  of  energy  equation,  better  known  as  Bernoulli's 

equation,  requires  that  for  any  two  points  1  and  2  ir  a  stream  tube: 

2  2 

_1_    +hi  +  zi=^?+h2  +  Z2  +  hL  (34) 

where 

V,  and  Vp  =  velocity  at  points  1  and  2,  respectively 
h.  and  h?  =  static  head  at  points  1  and  2,  respectively 
Z,  and  Z„  =  elevation  at  points  1  and  2,  respectively 
h,       =  the  total  specific  energy  lost  in  the  journey  of  the 
fluid  from  point  1  and  2 
Figure  1  shows  one  such  stream  tube.  The  velocity  at  point  1  is  very  low 
and  can  be  taken  to  be  zero.  The  elevation  is  the  same  for  both  points 
1  and  2  and,  thus,  cancel.  Loss  of  energy  as  water  moves  from  point  1 
to  2  was  neglected.  Also  since  the  static  heads  were  derived  from  the 
gage  pressures,  the  head  outside  the  tank  is  necessarily  zero.  With  these 
boundary  conditions,  Equation  34  becomes: 

V2  =  /~2gh  (35) 

V2  is  the  ideal  velocity  of  water  issuing  from  the  orifice  and  shall 
henceforth  be  referred  to  as  simply  V. 

Losses  due  to  various  factors  such  as  the  fluid  viscosity  make  the  in- 
troduction of  a  correction  factor,  Cw,  essential.  Hence;  V  =  Cv  /2gh  . 
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Figure  1.  Discharge  Through  an  Orifice  in  a  Tank. 


-^      Direction  of  Flow 


Figure  2.  Discharge  Through  an  Orifice  in  a  Pipe. 
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Continuity  requires  that: 

q'   =  aV  (36) 

where 

q'  =  orifice  discharge 

a  =  area  of  the  orifice 
Anywhere  outside  the  orifice,  the  true  area  of  discharge  is  always  less 
than  the  area  of  the  orifice.  Another  correction  factor,  Ca,  is  intro- 
duced to  take  care  of  the  reduced  area.  Thus,  the  true  discharge  becomes: 

q  =  CaCvA/TgTT  (37) 

The  coefficient  of  discharge,  C^,  is  then  defined  as  the  product  of  the 
two  correction  factors  (i.e.  Crf  =  C  L). 
Hence: 

q  =  CdA  /TgF  (38) 

Equation  38  has  been  used  widely  in  determining  orifice  flows  even 
in  the  case  of  manifold  pipes.  Nevertheless,  since  the  premises  upon 
which  the  theory  was  built  is  not  the  same  as  what  is  found  in  the  case 
of  a  manifold  pipe,  reassessment  of  the  various  factors  involved  is  in 
order. 

One  of  the  approximations  made  in  the  above  analysis  was  that  the 
velocity  at  point  number  1  was  zero  due  to  the  fact  that  the  amount  of 
water  that  exits  through  the  orifice  is  relatively  very   small  as  com- 
pared to  the  total  amount  of  water  in  the  tank.  In  the  case  of  a  manifold 
pipe,  however,  since  most  of  the  water  is  moving  past  the  orifice,  the 
velocity  inside  the  pipe  cannot  be  assumed  to  be  zero.  This  provides  a 
basis  for  another  new  approach  to  the  analysis. 
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If  Bernoulli's  equation  is  written  fcr  two  points,  one  inside  a  pipe 
directly  opposite  the  center  of  an  orifice  and  the  other  outside  on  an 
axis  of  discharge  through  the  orifice,   the  two  points  can  be  referred  to 
as  point  1  and  point  2,   respectively.     Bernoulli's  equation,  given  as 
Equation  34,  can  again  be  used. 

The  orifice  in  this  pipe  shall    in  this  analysis  be  taken  to  be  on 
the  side  of  the  pipe  such  that  the  axis  of  discharge  through  it  is  at 
right  angles  to  the  axis  of  flow  inside  the  pipe  and  the  discharge  exits 
horizontally.     The  elevation  of  the  two  points  1  and  2  are,   therefore, 
the  same.     Hence  the  terms  Z     and  Z~  in  Equation  34  are  eliminated.     In 
addition,  the  static  head  at  point  2  (i.e.   hp),   shall   also  be  zero  since 
the  gage  pressure  outside  the  pipe  is  zero. 

Equation  34,   then,   reduces   to: 


V2     =  /29   (  4?+  hl  "  hL}  ^9) 

If  the  head  loss  due  to  branching  and  exit  is  taken  to  be  a  fraction 
of  total  head  inside  the  pipe,  that  is  both  dynamic  and  static,  the  head 
loss  can  then  be  expressed  as; 

2 

2"g  ' 


Equation   39  then  becomes: 


V,2  V,2 


h-     7   29   C(4g     +  V    "   Khl?+hl^ 


(41) 
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Therefore: 


/         V  2 
V2  =     (1  -  K)h     Jig   (-L.    +  h2)  (42) 


If  we  let  Cy  =   (1  -  K) 2,  then: 


V2  =  Cy  >/2g   (  ^+  hj)  (43) 

If  the  total   area,  occupied  by  the  discharging  flow  while  inside 
the  pipe,   is  A',   then  continuity  requires  that: 

A'Vj  =   (CCA)V2  (44) 


where 


C     =  correction  factor  for  the  cross-sectional   area  of  the  discharge 
that  comes  out  of  the  orifice  since  a  contraction  always  occurs; 
the  factor  is  also  known  as  the  coefficient  of  contraction 

a  =  area  of  the  orifice 


Thus 

1 
Equation  43  becomes: 


C  A 
V,  =    4r    V,  (45) 


r  ccA  v22 

V2     =     Cy  V  2g     [(-Jr      -^    +  hj]  (46) 


Simplifying,  we  obtain: 


Cv     /2gh, 
\     ■  -  '    (47) 

7l  -  ( Vc  I  y- 

A  2 
The  factor,  (CyC  j,)  ,  in  the  demoninator  of  the  equation  is  so  small 

compared  to  1  that  it  can  be  assumed  to  be  zero.  In  which  case: 

V2  =  Cv  /2g7i7"  (48) 
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The  discharge  out  of  the  orifice  then  becomes: 


q  =  (CCA)  Cy  /  2ghx  (49) 

Thus: 

q  =  CdA  /TgtT  (50) 

where 

h  =  h1  in  the  preceding  equation 

A  completely  different  approach  to  the  analysis  of  orifice  discharge 
is  worth  considering.  In  the  interest  of  clarity,  the  following  condi- 
tions for  both  the  pipe  and  the  flow  through  it  shall  be  stated  as: 

1.  The  pipe  is  horizontal 

2.  The  orifices  are  in  the  vertical  plane 

3.  The  end  of  the  pipe  is  plugged 

With  these  conditions  we  can  then  move  on  to  say  that  at  any  point  of 
exit,  the  force  that  will  push  the  water  through  the  hole  causing  it  to 
change  directions  must  be  the  gage  pressure  inside  the  pipe.  The  angle 
that  the  discharge  out  of  the  orifice  makes  with  the  axis  of  flow  at  the 
point  of  exit  depends  upon  how  much  of  its  axial  momentum  it  has  lost. 
If  it  loses  all,  then  the  angle  has  to  be  90  .  However,  since  this  is 
seldom  the  case,  we  shall  assume  a  general  case  in  which  the  water  re- 
tains a  significant  amount  of  its  momentum  after  exit. 

Figure  2  shows  the  diagram  of  water  flowing  out  of  one  of  the 
orifices  of  a  pipe.  Since  the  water  is  coming  out  at  an  angle  6,  to  the 
axis  of  flow  inside  the  pipe,  the  actual  cross-sectional  area  of  the 
stream  shall  be: 

a  =  a  sin  6  (51) 

where 

a  =  the  area  of  the  orifice 
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It  can  further  be  stated  that  the  effective  force  acting  on  the  issuing 

stream  of  water  must  be  one  that  is  applied  in  the  same  direction  as  the 

flow.  Since  the  water  pressure,  P,  at  the  orifice  acts  in  all  directions, 

the  effective  force  shall  be  the  pressure  applied  on  the  cross-sectional 

area  of  the  water  stream.  Thus,  the  force,  F,  can  be  given  by: 

F  =  P  a  sin  9  (52) 

The  momentum  principle  requires  that: 

ZF  =  -j^r   material  pVdV  =  -r|  f  control  pVDV  +  [control  pV(V-ds)  (53) 
-V volume  J  volume        •'volume 

where  the  term  IF  on  the  left  represents  all  the  forces  acting  on  a  fluid 
while  the  terms  on  the  right  represent  the  total  rate  of  change  of  momen- 
tum. These  can  be  rewritten  as: 


if 


3t    JJ 


pV(V-ds) 


(54) 


where 
M 

9M 
3t 


VpdV,  the  momentum  of  fluid  within  the  control  volume 


for  the  fluid  initially  within  a  control  surface 


-  the  rate  of  change  of  momentum  within  the  control  volume 


pV(V-ds)  =  the  momentum  flow  rate  outward  through  the  control 


surface 


For  steady  flow,  that  is  when  the  velocity  at  any  one  point  does  not 


SM 


change  with  time,  the  |^  goes  to  zero  if  the  density  of  the  fluid  is  also 


assumed  to  be  constant.  Hence 


IF 


pV(V-ds) 


This  can  be  simplified  into  the  form: 

ZF  =  (ftV)  leaving  -  (MV)  entering 


(55) 
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Breaking  the  forces  into  their  various  components, 

ZF  =  (MVJ  leaving  -  (MVJ  entering  (56) 

X        X  A 

ZFV  =  (MVy)  leaving  -  (MV  )  entering  (57) 


EF2  =  (MV  )  leaving  -  (MV  )  entering  (58) 


where 


F  ,   F  ,   F     =  the  x,  y,  and  z  components  of  the  sum  of  the  forces 
x       y       z 

V   ,  V   ,  V     =  the  x,  y,  and  z  components  of  the  exit  and  inside 
velocities  of  the  issuing  water. 
In  Figure  2,  a  flow  was  shown  issuing  out  of  an  orifice  at  an  angle 
6  to  the  axis  of  flow.     If  the  effective  force  acting  on  it  is  P  a  sine 
as  earlier  stated,  then,   using  the  convention  indicated  in  the  figure, 
the  y-component  of  the  force,   F   ,  d,    is  given  by: 

F     =  Pasin26  (59) 

Since     F       is  the  only  force  acting  in  the  y  direction  and  the  initial 
velocity  of  the  water  in  the  same  direction  was  zero,   the  principle  of 
momentum  requires  that: 


Fy  -   pq    (Ve  sine   -  0)  (60) 


where 


V     =  the  exit  velocity 
e  J 

q     =  the  discharge  out  of  the  orifice 
p    =  water  density 

Thus: 

2 
Pasin  6  =  qV  sine 

But,   since  asin6  is  the  actual   cross-sectional   area  of  the  issuing 

stream: 

ve=    «?R5  <61» 
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n  0 

'"•   pq  '   asine     '   sinQ  =  Pasin  6 

2  2       2 

or  pq     =  Pa  sin  6 

But  P  =  pgh 

where 

h  =  head  in  feet  of  the  pressure  at  the  orifice 

2  2       2 

. .   pq     =  gha  sin  6 

2  2       2 

or     q     =  gha  sin  6 

.-.     q     =    SM    .     a  ^  (62) 

/  2 

This  equation  is  of  the  same  form  as  Equation  38.  So,  if  we  are  to 

use  the  latter,  then  we  shall  know  that  C,  =  ^2M.   , 

d  fY 

The  above  analysis  indicates  that  the  discharge  coefficient  is  a 
function  of  the  angle  of  discharge  for  the  orifice.  It  can  be  recalled 
that  the  angle  of  discharge  itself  has  been  stated  to  be  a  function  of 
both  the  velocity  of  flow  inside  the  pipe  and  the  pressure  at  the  orifice. 
It  can  therefore  be  inferred  that  the  discharge  coefficient,  C  ,,  is  a 
function  of  both  the  velocity  of  the  fluid  as  it  passes  the  orifice  and 
the  pressure  head  at  the  same  point. 

If  there  is  no  analytical  way  of  relating  the  angle  of  discharge  to 
the  two  factors  affecting  it,  then  an  empirical  relationship  can  be 
sought  through  experiments  to  determine  a  direct  relationship  between 
the  discharge  coefficient  and  the  two  factors. 


31 
Materials  and  Procedure 

The  experiment  was  conducted  in  the  hydraulics  laboratory  of  the 
College  of  Engineering  at  Kansas  State  University.  Some  of  the  basic 
equipment  used  was  already  available.  An  underground  tank  in  the  labora- 
tory was  employed  as  the  water  source. 

Among  the  available  equipment  was  a  horizontal  centrifugal  pump 
driven  by  a  20  h.p.  variable  speed  motor  shown  in  Figure  3.  This  was 
used  for  drawing  the  water  from  the  tank.  The  pump  was  primed  by  means 
of  a  vacuum  pump  which  can  be  seen  in  Figure  4. 

The  equipment  was  initially  arranged  such  that  water  was  pumped  into 
a  six-inch  pipe  which  led  to  a  six-inch  valve.  The  valve,  shown  in 
Figure  5,  controlled  the  rate  of  flow.  Flow  then  passed  through  an- 
other short  section  of  six-inch  pipe  into  a  six-inch  to  eight-inch  pipe 
adapter.  A  90  elbow  was  connected  to  the  adapter  and  to  this  elbow 
was,  in  turn,  connected  a  tee.  The  middle  outlet  of  the  tee  was  plugged 
so  that  the  flow  that  came  out  of  the  elbow  could  pass  straight  through 
it.  It  was,  thus,  used  simply  as  a  short  section  of  pipe.  From  the  tee, 
the  flow  passed  into  a  second  90°  elbow.  Another  connection  of  20-foot 
length  of  an  eight-inch  plain  pipe  was  made  at  the  other  end  of  the 
elbow.  To  this  pipe  was  connected  the  test  pipe  as  shown  in  Figure  6. 

Later,  when  surging  was  found  to  be  a  problem,  the  set  up  was  modi- 
fied by  making  the  flow  pass  through  a  stand  pipe.  The  tee  used  earlier 
as  a  pipe  section  was  connected  directly  to  the  six-to-eight-inch  adapter 
mentioned  earlier.   The  other  end  of  it  was  plugged  while  a  plain  pipe 
about  10  feet  in  length  was  connected  to  its  middle  outlet  and  the  flow 
was  made  to  go  vertically  upwards  through  the  pipe.  The  flow  then  dis- 
charged into  the  stand  pipe  after  passing  through  two  90°  elbows  which 
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Figure  3.     Pump  in  Kansas  State  University 
Hydraulics  Laboratory. 


Figure  4.  Vacuum  Pump  for  Priming  Pump  and 
Associated  Valves  and  Piping. 
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Figure  5.  The  Flow  Control  Valve. 


Figure  6.  Adjustable  Wooden  Supports  that  Held  the  Pipe. 
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were  connected  together.  The  20-foot  length  of  pipe  mentioned  earlier 
was  then  directly  connected  to  the  stand  pipe. 

The  test  pipe  used  was  thirty  feet  long  and  eight-inches  in  diameter. 
The  eight-inch,  20-foot  long  pipe,  was  used  for  stabilizing  the  flow  as 
it  came  out  of  the  last  elbow  before  it  entered  the  test  section.  Measure- 
ments of  the  total  discharge  passing  through  the  pipe  were  made  with  a 
Jaccuzi  flow  meter.  Figure  7  shows  the  flow  meter.  A  differential  mano- 
meter, shown  in  Figure  8,  was  used  whenever  the  determination  of  head 
difference  between  two  points  was  desired.  Copper  tubings,  about  1  1/2 
inch  in  length  and  0.125-inch  internal  diameter,  were  glued  to  the  test 
pipe  at  the  various  orifices  and  through  them  holes  were  drilled  into 
the  pipe  to  provide  pressure  tappings  as  shown  in  Figure  9.  Head  read- 
ings were  taken  from  a  set  of  manometers  mounted  in  a  protective  steel 
frame  against  a  set  of  scales  marked  from  zero  to  seventy-two  inches. 
These  are  shown  in  Figure  10.  Connections  between  the  manometers  and 
the  copper  tubings  were  made  with  0.25-inch  inside  diameter  flexible 
polyethelene  tubings. 

An  initial  test  was  made  to  determine  the  friction  coefficient  of 
the  test  pipe.  Using  two  pressure  tappings,  27  1/2  feet  apart,  with 
one  being  15-inches  from  the  inlet  end  and  the  other  being  the  same 
distance  from  the  outlet  end,  the  head  loss  between  the  two  points  was 
determined  with  the  aid  of  the  differential  manometer.  The  velocity  of 
flow  was  measured  with  the  flow  meter.  Readings  were  taken  at  ten  dif- 
ferent points  across  a  section  of  the  pipe,  each  representing  one  of 
ten  equal  areas  bounded  by  ten  concentric  circles.  The  average  of  the 
ten  was  then  taken  to  be  the  average  velocity. 


35 


,i-:C.-> 


i. 

O! 

+J 

E= 
o 

c 

10 


to 


s- 

QJ 


cc 

i- 
3 
cn 


W 
+J 

J! 

O 


N 
N 

(J 

u 

■"3 

0) 


3 
C31 


36 


Figure  9.  Pressure  Tappings  on  the  Test  Pipe. 


Figure  10.  Set  of  Manometers 
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The  test  pipe  was  then  disconnected  from  the  20-foot  pipe.  Holes 
were  drilled  in  the  pipe  with  a  drill  press.  During  the  drilling  pro- 
cess, the  pipe  was  placed  in  a  wooden  saddle  such  that  about  half  of  its 
perimeter  was  in  the  saddle.  A  metal  band  was  strapped  around  the  upper 
part  of  the  pipe  and  screwed  onto  the  saddle  for  added  stability.  The 
drill  bit  was  lowered  gently  onto  the  pipe,  each  time,  and  was  run  at 
its  slowest  speed  to  avoid  any  wobbling  that  would  easily  jeopardize 
the  preciseness  of  the  holes.  The  arrangement  can  be  seen  in  Figure  11. 
Only  six  holes  were  drilled  on  the  test  pipe,  each  13/16  inch  in  diameter 
and  all  equally  spaced  at  5  feet  from  center  to  center.  A  manually 
operated  crane  was  used  for  lifting  the  barrels  which  were  used  for  col- 
lecting the  discharge  out  of  the  various  orifices  of  the  pipe.    The 
crane  used  is  shown  in  Figure  12. 

The  pipe  was  then  connected  back  onto  the  20-foot  pipe  and  flow 
diverting  devices  were  mounted  on  it  with  one  unit  at  each  hole.  The 
open  end  was  plugged,  thus  becoming  a  dead  end. 

Water  was  run  through  the  system  and  discharged  through  the  six 
orifices.  The  static  head  at  each  outlet  with  reference  to  the  center 
line  of  the  orifice  was  recorded  from  the  set  of  manometers  previously 
mentioned.  The  discharged  water  was  collected  in  barrels  for  some  re- 
corded time  interval.  The  barrels  were  always  weighed  while  empty,  just 
before  the  water  collection  started.  They  were  reweighed  again  while 
containing  the  water. 

The  flow  diverting  device  mounted  on  the  pipe  initially  was  one 
that  allowed  the  discharge  out  of  all  the  orifices  to  be  diverted  simul- 
taneously. The  device  consisted  of  six  identical  units,  all  connected 
together  and  activated  from  the  same  point  at  one  time.  Each  unit  was 
made  up  of  three  pieces  of  metal  bars  as  illustrated  in  Figure  13.  One 
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Figure  11.  Arrangement  for  Drilling  the  Orifices 


SWWSWP 


Figure  12.     Manually  Operated  Crane  for 
Lifting  a  Can-full   of  Water. 


39 


O 


s 


10' 


3" 


Figure  13.  Components  of  the  Flow  Diverting  Device: 


Figure  14.  The  Flow  Diverting  Device  in  Place, 
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about  10  inches  in  length  while  the  other  two  were  7  inches  and  5  inches. 
In  the  interest  of  clarity,  the  three  pieces  will  be  referred  to  as  A, 
B,  and  C.  The  pieces  B  and  C  were  welded  together  forming  a  right  angle 
at  their  joint.  On  the  piece  C,  a  hole  was  drilled  about  4.5  inches  from 
the  joint.  Another  hole  was  drilled  on  piece  A,  1/2  inch  from  one  end. 
Through  these  holes  the  two  pieces,  A  and  C,  were  connected  by  means  of 
a  bolt  which  acted  as  a  pivot  point  around  which  piece  A  could  swing  freely. 
Two  more  holes  were  drilled  in  piece  A,  one  at  each  end.  A  metal  plate 
was  bolted  onto  it  at  the  end  furthest  from  the  pivot  point  while  two 
handles  were  connected  to  the  other  end  such  that  a  pull  on  one  of  them 
could  cause  the  piece  to  swing  one  way  or  the  other.  Thus,  when  the 
device  was  put  in  place  around  an  orifice,  the  plate  could  be  made  to 
either  cut  across  the  flow  or  be  swung  away  from  it.  The  device  was 
welded  to  a  pipe  clamp  and  then  mounted  on  the  pipe  around  an  orifice. 
Figure  14  shows  the  flow  diverting  device  in  place.  All  the  handles 
were  connected  together  such  that  any  action  at  one  end  could  cause  a 
simultaneous  identical  response  in  all  units. 

This  was  attained  by  turning  the  pipe  so  that  the  discharge  exited 
horizontally.  When  the  barrels  were  ready  to  be  put  in  place,  the  diver- 
ter  was  activated  so  that  the  plates  at  each  unit  cut  across  the  flow 
and  remained  there,  thus  obstructing  flow.  The  barrels  were  then  placed 
in  position  and  the  diverter  was  again  activated,  this  time  to  remove 
the  plate  obstruction.  At  the  same  time,  a  stop  watch  was  started.  When 
the  water  in  the  barrel  collected  up  to  a  selected  level,  the  diverter 
was  activated  for  a  third  time  to  cut  the  flow  into  the  barrels.  Weigh- 
ing then  followed  as  mentioned  earlier.  The  process  was  repeated  for 
four  different  flow- rates.  Thus  four  runs  were  made. 
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The  test  pipe  was  then  removed  for  a  second  time.  The  holes  were 
enlarged  to  1  1/4  inches  in  diameter  by  using  a  drill  bit  of  proper  size. 
The  procedure  used  for  the  drilling  was  very   much  the  same  as  that  de- 
scribed earlier  in  the  case  of  the  smaller  holes.  Additional  care  was, 
however,  taken  to  make  sure  that  the  bit  was  well  centered. 

The  pipe  was  then  put  back  in  place.  Water  was  run  through  as  before 
only  this  time,  the  pipe  was  oriented  such  that  the  orifices  discharged 
directly  downwards  to  facilitate  collection  of  the  flow.  A  different 
diverter  had  to  be  used  because  the  first  diverter  could  not  control 
flow  through  the  enlarged  orifices.  The  second  diverter  consisted  of 
a  single  unit  and  was  made  up  of  two  pieces  of  rigid  plastic  pipes  con- 
nected together  with  a  75°  elbow.  A  chain  tied  around  the  shorter  of 
the  two  pieces  was  made  to  go  around  the  test  pipe  at  the  required  posi- 
tion and  the  loose  end  was  hung  on  to  a  hook  fastened  to  the  part  of  the 
chain  that  was  tied  to  the  plastic  pipe.  The  arrangement,  which  is  shown 
in  Figures  15  and  16,  was  used  to  divert  water  into  a  barrel  for  a  meas- 
ured time  interval.  This  interval  was  recorded  with  an  electronic  stop- 
watch. So,  at  each  orifice  the  diverter  was  swung  into  place  and,  simul- 
taneously, the  start  button  of  the  stop  watch  was  pressed.  The  water 
flowed  through  the  plastic  tubing,  into  the  elbow  and  through  the  second 
tubing  into  the  barrel  which  was  sitting  on  the  scale.  As  in  the  case 
of  the  first  set  of  orifices,  in  this  also,  four  runs  were  made  by  vary- 
ing the  flow  rate  into  the  system. 

The  diverter  was  swung  away  after  some  recorded  time  interval  and 
readings  were  taken  on  the  scale.  The  same  procedure  was  followed  for 
determining  the  discharge  out  of  the  rest  of  the  orifices  taking  one 
outlet  at  a  time.  Four  runs  were  made  with  this  arrangement. 
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The  scale  used  was  a  regular  platform  scale  that  read  to  one-tenth 
of  a  pound.  However,  an  estimate  to  the  nearest  hundredth  was  acceptably 
possible.  Thus  weight,  instead  of  volume,  was  measured  for  each  sampled 
discharge.  Both  volume  and  the  average  velocity  at  various  sections 
could  be  determined. 
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Results 

The  test  for  the  determination  of  the  friction  coefficient  for  the 
eight-inch  diameter  test  pipe  showed  that,  when  the  average  velocity 
inside  the  pipe  was  4.255  feet  per  second,  the  total  head  loss  between 
the  two  tappings,  27  1/2  feet  apart,  was  2.8  inches.  The  friction  coef- 
ficient, Cm,  was  determined  using  the  Hazen-Williams '  equation  of  velocity. 
The  equation  is: 

V  =  1.318CHR°-63S0-54  <63> 

where 

V  =  velocity  in  the  pipe,  cfs 
Cn  =  pipe  friction  coefficient 
R  =  the  pipe  hydraulic  radius,  ft 
S  =  the  hydraulic  slope 
The  head  loss  due  to  fricition  can  be  expressed  as: 

hf  =  SL  (64) 

where 

h-  =  the  head  loss  due  to  friction,  ft 

L  =  length  of  segment  of  the  pipe  under  consideration;  ft 
After  solving  Equation  64  for  S,  substituting  in  Equation  63  and  rearrang- 
ing, the  following  expression  can  be  established: 

1 


CH  ■       [       V     Q  fi~]       A-]   1,8S  (65) 

H  1.318R0"63  hf 

Using  Equation  65  the  friction  coefficient  was  calculated  to  be  131.48. 

Following  the  determination  of  the  friction  coefficient,  two  sets  of 
tests  were  carried  out  using  the  same  test  pipe.     These  tests  were  con- 
ducted after  drilling  holes  in  the  pipe.     The  aim  as  stated  in  the  objec- 
tives was  to  determine  the  flow  rate  from  the  orifices  and  the  static  head 
at  each.     The  approach  velocity  toward  the  orifices  and  the  discharge 
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coefficient  at  each  were  then  calculated.   Holes  of  13/16-inch  diameter 
were  used  in  the  first  set  of  tests  while  1.25-inch  diameter  holes  were 
used  in  the  second  set.  The  results  are  shown  in  Tables  1  and  2,  respec- 
tively. 

Suggestions,  from  earlier  analyses  and  experiments,  indicated  that 
the  coefficient  of  discharge  is  affected  by  both  the  velocity  approach 
and  the  head  at  an  orifice. 

Regression  analyses  using  a  computer  were,  therefore,  carried  out  to 
determine  the  unknown  effects  and  the  full  dimension  of  the  nature  of 
any  relationship  that  may  exist  between  them. 

The  computer  program  used,  known  as  the  "Statistical  Analysis  System," 
or  simply  SAS  performs  various  statistical  analyses  depending  upon  the 
procedure  required  and  the  model  indicated. 

Statistical  theory  suggested  that  the  relationship  between  any  two 
or  more  quantities  could  be  expressed  as  a  polynomial  equation.  So  if 
a  set  of  experimental  values  of  two  or  more  quantities  were  available, 
in  order  to  determine  if  some  relationship  existed  between  them,  attempts 
could  be  made  to  express  one  of  the  quantities  as  a  polynomial  function 
of  the  rest.  However,  a  polynomial  could  be  of  the  first,  second,  third, 
or  any  other  order.  But  in  choosing  an  order  for  a  model,  it  is  always 
logical  to  choose  the  first  order  model  as  the  starting  one. 

For  the  analysis  of  the  data  obtained  in  this  study,  the  regression 
procedure  of  the  SAS  program  was  used.  The  first  order  model  was  tried 
first,  then  the  second,  then  the  third  and  so  on.  For  every   order  of 
model  considered  an  F-test  was  conducted  to  determine  if  the  order  of 
the  model  in  question  had  shown  a  significant  reduction  of  the  residual 
sum  of  the  squares,  at  5%   level,  as  compared  to  the  model  order  immedi- 
ately preceding  it.  For  instance,  an  F-test  was  made  to  find  out  if  the 
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second  order  model  had  shown  a  significant  reduction  of  the  residual  sum 
of  the  squares  over  the  first  model.  When  two  consecutive  models  show 
insignificant  reductions  of  sum  of  squares  over  their  preceding  ones, 
then  the  model  immediately  preceding  the  two  was  taken  as  the  model  of 
best  fit.  For  example,  if  the  second  order  model  showed  no  significant 
reduction  of  the  residual  sum  of  the  squares  over  the  first,  and  nor  did 
the  third  over  the  second,  then  according  to  this  methodology,  the  first 
order  model  was  taken  as  the  model  with  the  best  fit. 

Using  the  above  methodology  in  the  analysis  in  this  study,  a  fourth 
order  model  was  obtained  having  15  terms.  Since,  some  of  the  terms  were 
of  no  significance  in  the  relationship  between  the  quantities,  another 
procedure  in  the  SAS  program  was  used  to  determine  the  significant  terms 
in  the  relationship.  The  procedure,  known  as  the  "stepwise  procedure," 
took  only  one  term  initially  and  added  the  other  terms  in  the  model  one 
at  a  time.  An  F-test  was  made  in  each  case  to  determine  the  significance 
of  the  term  to  the  model.  Thus,  all  insignificant  terms  were  rejected. 
Results  showed  that  the  discharge  coefficient  at  any  orifice  could 
be  expressed  as  a  third  order  model  and  a  function  of  both  the  velocity 
of  approach  and  the  pressure  head  at  the  orifice.  The  prediction  equation 
obtained  was: 

Cd  =  0.5836  +  0.3723V  -  0.01098IW  -  0.346V2  +  0.1084V3       (66) 
where 

C  .  =  discharge  coefficient  of  an  orifice 

V  =  approach  velocity  toward  the  orifice 

h  =  pressure  head  at  the  orifice 

The  appearance  of  the  pressure  head  in  only  one  term  and  the  pre- 
dominance of  the  approach  velocity  term  in  the  equation  raised  some 
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questions  about  the  significance  of  the  pressure  head  effect  upon  the 
coefficient. 

A  similar  analysis  was  carried  out,  taking  the  approach  velocity  as 
the  only  factor  of  concern.  Again  another  prediction  equation  of  the 
same  third  order  was  obtained.  It  was  established  as: 

Cd  =  0.5883  +  0.3106V  -  0.3141V2  +  0.0898V3  (67) 

A  comparison  of  discharge  coefficients  predicted  by  each  of  these 
two  equations  with  the  observed  values  showed  a  good  agreement.  The 
comparisons  are  given  in  Tables  3  and  4.  The  differences  between  the 
observed  and  the  predicted  values  were  found  to  be  statistically  insig- 
nificant, after  performing  an  F-test.  It  can  be  recalled  that  part  of 
the  objective  of  this  study  was  to  determine  if  the  Enger  and  Levy 
equation,  given  earlier  as  Equation  3,  can  be  applied  to  the  case  of 
eight-inch  irrigation  pipes.  It  was  found  that  the  accuracy  of  its  pre- 
diction depended  very  much  on  how  good  an  assumption  was  made  of  the  dis- 
charge coefficient  for  the  dead  end.  Values  predicted  by  this  equation 
can  also  be  seen  in  the  comparison  Tables  3  and  4. 
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Discussion 

The  results  of  the  theoretical  discussion  and  analysis  in  this  study 
suggested  that  the  discharge  coefficient  for  any  orifice  is  a  function 
of  the  angle  at  which  the  discharge  exits  from  the  orifice  and  was  given 
in  Equation  62.  It  was  further  suggested  that  the  angle  of  discharge 
itself  is  a  function  of  both  the  velocity  of  approach  towards  and  the 
pressure  at  the  orifice.  The  obvious  inference  was  therefore  that  the 
discharge  coefficient  is  a  function  of  the  approach  velocity  and  the 
static  head  at  the  orifice.  The  analysis,  however,  failed  to  establish 
any  relationship  similar  to  those  obtained  as  a  result  of  this  experi- 
ment (Equations  66  and  67).  These  equations,  though  subject  to  some 
limitations,  had  nonetheless  strengthened  the  claim  that  was  based  on 
the  theoretical  analysis.  The  limitation  alluded  to  here  is  one  of 
scope  of  validity. 

It  must  be  conceded  that  this  study  was  not  designed  to  include 
many  other  factors  that  are  likely  to  play  a  role  in  the  influence  of 
the  discharge  coefficient.  These  factors  were  dropped  from  consideration 
on  the  grounds  that  some  earlier  studies  had  found  their  effect  insigni- 
ficant as  compared  to  the  effects  of  approach  velocity  and  static  head. 
Nevertheless,  since  a  generalization  cannot  be  made  without  some  risk 
of  oversimplification  of  the  problem,  the  decision  to  consider  those 
factors  ineffective  is,  to  some  extent,  tantamount  to  an  assumption. 
This  exposes  one  of  the  sources  of  limitations  of  the  obtained  empirical 
relations. 

Moreover,  the  prediction  Equations  66  and  67  may  be  valid  only 
within  the  range  of  the  data  collected.  Even  though  this  can  best  be 
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ascertained  by  a  comparison  between  values  predicted  by  the  equation,  and 
those  calculated  from  some  observed  data  of  flow  well  beyond  the  range 
recorded  in  this  experiment. 

This  notwithstanding,  the  empirical  formula  can  be  used  to  determine 
the  rate  of  flow  out  of  the  various  orifices  of  a  pipe  when  the  total 
discharge  and  static  head  at  the  point  of  entry  are  known.  The  accuracy 
of  the  values  will  of  course  depend  on  how  closely  the  conditions  of  the 
system  in  question  resemble  those  under  which  this  experiment  was  conducted. 
Enger  and  Levy's  Equation  showed  good  agreements  with  the  observed  data 
in  some  of  the  runs.  However,  since  the  value  of  any  coefficient  can 
be  predicted  only  after  assuming  a  value  for  the  coefficient  of  the  ori- 
fice nearest  the  dead  end,  the  observed  value  for  it  was  used.  Hence  the 
value  of  the  observed  and  that  obtained  by  the  prediction  equation  are 
always  the  same. 

In  order  to  find  out  how  good  a  prediction  for  orifice  discharge  can 
be  made  employing  those  empirical  relations,  the  data  of  a  different  ex- 
periment was  used  for  comparison.  The  experiment  was  performed  by  Spomer 
(1969)  and  the  data  gave  the  total  discharge  into  a  manifold  pipe  along 
with  the  static  head  at  the  entry  point.  Also  included  were  the  observed 
discharge  out  of  the  various  orifices. 

The  following  methodology  was  used  for  calculating  the  orifice  dis- 
charges using  the  empirical  equations  66  and  67. 

First,  the  inflow  and  the  head  at  entrance  to  the  gated  pipe  was 
known.  The  velocity  of  approach  towards  the  first  orifice  from  the  inlet 
end  was  then  determined  using  the  continuity  equation: 

Q  =  AV  (68) 


5S 


where 

Q  =  the  total   flow  at  entry  into  the  pipe 

A  =  the  cross-sectional  area  of  the  pipe 

V  =  the  velocity  inside  the  pipe 
Then,  using  the  empirical  equation  66,  the  discharge  coefficient  for  the 
first  orifice  was  obtained  under  the  assumption  that  the  pressure  head 
just  before  the  orifice  is  the  same  as  that  at  the  orifice.     This  pres- 
sure head  was  determined  by  subtracting  the  head  loss  due  to  friction 
between  the  point  of  measurement  of  head  at  entry  and  the  orifice. 
Hazen-Williams '  equation  of  velocity  was  used  in  this  connection. 

Since  the  head  at  the  orifice  was  known  along  with  the  discharge 
coefficient,  the  flow  out  of  the  first  orifice  was  calculated  using 
Equation  38.       Bernoulli's  equation  was,  therefore,  written  between  the 
point    just  prior  to  the  orifice  and  one  immediately  following  it.     Thus: 

V  2  V2 

hl  +  2i~    *     h2  +  W  m 

where 

h,  =  the  static  head  immediately  before  the  orifice 
V,  =  approach  velocity  toward  the  orifice 
h„  =  the  static  head  immediately  after  the  orifice 
V2  =  velocity  of  water  leaving  the  first  orifice 
The  velocity  of  water  inside  the  pipe  leaving  the  orifice,  V2,  can 
be  expressed  as: 

V  Vj  -  Ji  (70) 

This  equation  is  true,  since  the  total  flow  inside  the  pipe  after  passing 
the  first  orifice,  Q2,  is  equal  to  Q.  -  q,.  Another  form  of  this  state- 
ment from  which  the  above  equation  was  directly  derived  can  be  expressed 
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as,  AV?  -  AV.   -  q. .     The  energy  equation,   therefore,  became: 

ql     2 

v 2     <vi  -  r  > 

h2  =  h1+^ (71) 

From  this  equation  the  static  head  immediately  after  the  orifice, 
h2>  was  calculated  and  the  velocity  V?  in  the  section  downstream  from 
the  first  orifice  was  calculated  from  Equation  70  .  This  velocity  became 
the  approach  velocity  toward  the  second  orifice.  The  static  head  immedi- 
ately before  the  second  orifice,  h^,  was  the  static  head  just  after  the 
first  orifice  less  the  friction  loss  in  between.  Again  Hazen-Williams ' 
equation  was  used  for  determining  the  head  loss  due  to  friction. 
Therefore: 

h2  =  hl  "  hfl-2  (72) 

where 

hi  =  static  head  just  before  the  second  orifice 

hfl-2  =  the  head  ''oss  due  t0  'frrict'ion  between  orifice  1  and  2 
Then,  onward,  the  process  was  repeated  until  the  dead  end  was  reached. 

Approximations  made  for  the  static  head  at  an  orifice,  in  this  pro- 
cedure, contain  the  assumption  that  the  static  head  just  before  the 
orifice  was  the  same  as  the  head  at  the  orifice.   Since  this  is  anything 
but  the  actual  case,  adjustments  were  made  to  correct  for  the  assumption. 
With  all  heads  both  before  and  after  each  orifice  calculated,  the  average 
of  the  two  was  taken  to  be  the  actual  head  at  the  orifices.  However, 
this  adjustment  was  started  with  the  first  orifice  from  the  inlet  end 
while  the  rest  were  left  as  they  were.  The  discharge  out  of  the  first 
orifice  was  recalculated  using  the  new  value  for  the  head.  The  process 
of  recalculating  other  heads  both  upstream  and  downstream  of  the  orifices 
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followed  in  accordance  with  the  methodology  described  earlier.  The  new 
value  for  the  head  was  examined  to  determine  the  percentage  difference 
between  it  and  the  previous  value  of  the  head.  When  the  per- 

centage was  found  to  be  less  than  1,  the  adjustment  process  proceeded 
to  the  next  orifice.  The  procedure  was  repeated  until  the  last  orifice 
was  reached.  Thus,  the  head  at  every  orifice  and  the  discharge  from 
each  were  established.  Table  5  compared  the  value  calculated  and  those 
actually  observed  by  Spomer.  Even  though  satisfactory  values  of  head 
everywhere  along  the  pipe  were  obtained,  the  computed  discharge  out  of 
the  orifices  were  not  equal.  In  order  to  attain  uniform  flow,  an  amount 
of  discharge  equal  to  the  quotient  of  the  total  inflow  into  the  pipe 
divided  by  the  total  number  of  the  orifices  had  to  issue  from  each  orifice. 
This,  then,  was  regarded  as  the  required  flow.  A  head,  other  than  the  one 
computed,  was  needed  to  produce  this  required  flow  out  of  each  of  the 
orifices.  When  the  computed  head  was  greater  than  the  required  one,  the 
difference  was  regarded  as  positive  and  when  it  was  the  other  way  around, 
the  difference  was  considered  negative.  All  the  differences  were  deter- 
mined associating  the  appropriate  signs  with  each. 

The  pipe  under  discussion  was  level  and  all  the  orifices  in  it  were 
at  the  same  elevation.  When  the  head  differences  were  subtracted  from 
the  pipe  elevation  at  e\/ery   orifice,  a  different  elevation  resulted  in 
each  case.  However,  since  the  pipe  could  not  be  bent  all  over  such  as 
to  obtain  the  ideal  elevation  at  each  orifice,  it  was  found  most  appro- 
priate to  determine  an  optimum  pipe  slope  that  would  approximate  most 
closely  the  ideal  elevations.  The  slope  of  the  line  of  least  squares  was 
calculated  and  was  taken  to  be  the  optimum  slope.  A  computer  program 
was  written  and  used  for  all  the  calculations  described  in  this  procedure. 
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Table  5.       Comparison  of  Calculated  and  Observed 
Discharge  Using  Spomer's  Data. 


Calculated 
Discharge 

Observed 
Discharge 

Difference 
of  the  Two 

Percentage* 
Deviation 

5.09 

5.00 

-0.09 

-0.02 

5.10 

4.91 

-0.19 

-0.04 

5.13 

5.02 

-0.11 

-0.02 

5.16 

5.06 

-0.10 

-0.02 

5.18 

5.04 

-0.14 

-0.03 

5.19 

5.05 

-0.14 

-0.03 

5.20 

5.01 

-0.19 

-0.04 

5.20 

5.02 

-0.18 

-0.04 

5.20 

5.06 

-0.14 

-0.03 

5.18 

5.00 

-0.18 

-0.04 

5.15 

5.04 

-0.11 

-0.02 

5.11 

5.01 

-0.10 

-0.02 

5.06 

5.01 

-0.05 

-0.01 

5.00 

5.09 

0.09 

0.02 

4.93 

5.20 

0.27 

0.05 

4.84 

5.24 

0.40 

0.08 

4.74 

4.99 

0.25 

0.05 

4.64 

5.02 

0.38 

0.08 

Mean  =  -0.02   Square  of  Standard  Deviation  =  0.04   T  =  -0.40 
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A  printout  obtained  as  a  result  is  shown  in  the  Appendix.  It  should  be 
pointed  out  here  that  the  program  was  written  such  as  to  allow  the  use 
of  it  by  anybody,  in  any  design  for  optimum  pipe  slope  that  would  provide 
uniform  flow,  so  long  as  the  total  inflow  and  the  head  just  before  the 
first  orifice  were  known. 

The  Enger  and  Levy's  equation  was  found  unworkable  under  the  method- 
ology just  described. 
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CONCLUSIONS 

The  following  were  concluded  as  a  result  of  this  study: 

1.  The  coefficient  of  discharge  of  an  orifice  in  the  side  of  a 

multiple  outlet  pipe  can  be  expressed  as  a  function  of  either  the 

approach  velocity  towards  the  orifice  or  as  a  function  of  both  the 

approach  velocity  and  the  static  head  at  the  orifice.     Further,  the 

value  of  the  coefficient  can  be  predicted  by  either  of  the  following 

two  equations: 

C,  =  0.5836  +  0.3723V  -  0.01098hV   -  0.346V2  +  0.1084V3 
d 

and 
Cd  =  0.5883  +  0.3106V  -  0.3141V2  +  0.0898V3 

2.  The  Enger  and  Levy  Equation  can  be  applied  to  the  case  of  a 
multi-outlet  eight-inch  irrigation  pipe  with  the  accuracy  of  pre- 
diction depending  to  a  large  extent  on  the  accuracy  of  the  assumption 
for  the  value  of  the  discharge  coefficient  of  the  orifice  at  the 
dead  end. 

3.  The  prediction  equations,  obtained  as  a  result  of  the  investi- 
gation, predicted  satisfactory  values  for  orifice  discharge  in  the 
case  of  different  pipe  and  orifice  sizes.     The  predicted  values,  also, 
compared  favorably  with  observed  values  obtained  from  a  different 
experiment. 

4.  A  methodology  was  established  for  determining  an  optimum  slope 
needed  to  provide  uniform  flow  out  of  various  orifices  of  a  gated 
pipe.     A  computer  program  was  developed  for  easy  application.     In 
reality  exact  uniformity  may  not  be  achieved  but  deviations  therefore 
should  be  practically  insignificant.     Lastly,  the  Enger  and  Levy 
Equation  cannot  work  with  this  methodology. 
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SUMMARY 

Non- uniformity  of  flow  through  the  various  orifices  of  a  gated  pipe 
irrigation  system  has  presented  a  problem  that  is  becoming  of  increasing 
concern.  At  present  farmers  try  to  overcome  it  by  adjusting  the  individual 
gates  until  equal  flow  is  approximated  as  nearly  as  possible.  Increasing 
costs  of  operation  due  to  high  labor  inputs  in  many  irrigation  systems 
call  for  automation  of  the  system.  If,  however,  gates  have  to  be  ad- 
justed manually  employing  both  the  skill  and  judgement  of  the  irrigator, 
any  automation  shall  be  far  from  complete. 

This  study,  therefore,  addressed  itself  to  the  problem.  A  review  of 
the  literature  carried  out  showed  that  other  people  have  looked  into 
similar  problems  in  studies  with  gas  and  air  systems.  Yet  others  have 
worked  with  multiple  outlet  pipes,  some  with  perforated  pipes  and  some 
even  looked  directly  at  the  problem  of  uniformity  of  flow  from  gated  ir- 
rigation pipes.  There  appeared,  generally,  to  be  an  agreement  between 
the  people  that  considered  the  problem  that  the  discharge  coefficient 
at  each  orifice  is  affected  by  the  approach  velocity  and  the  head  at  each 
orifice.  The  coefficient  was  accepted  as  a  factor  having  a  direct  in- 
fluence on  the  discharge  out  of  the  orifices.  Some  have  alluded  to  the 
possibility  of  the  effects  of  such  factors  as  the  orifice  to  pipe  diameter 
ratios  and  the  curvature  of  the  pipe  just  to  mention  a  few.  Others  have 
dismissed  them  as  insignificant  compared  to  the  effects  of  the  two  factors 
previously  mentioned.  Despite  all  these  suggestions,  none  but  one  re- 
lationship between  the  discharge  coefficient  and  the  approach  velocity 
and,  or  the  static  head  was  found  to  have  been  established.  Enger  and 
Levy  gave  that  relationship  expressing  the  discharge  coefficient  at  an 
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orifice  as  a  function  of  the  dynamic  and  static  head  upstream  of  an  orifice 
along  with  discharge  coefficient  at  the  dead  end.     This  latter  quantity 
has  to  be  assumed  every  time  the  equation  is  used.     The  equation  was  ob- 
tained in  a  study  with  a  2-inch  pipe  having  3/8-inch  orifices  and  was 
expressed  as:  2 

r     =   (i  _  -JS — )   r 
Ld       u  h     ;  Ld 

o 

An  investigation  was  carried  out  with  the  following  objectives: 

1.  To  determine  a  direct  relationship  between  the  orifice  dis- 
charge coefficient  and  the  parameters;  approach  velocity  and  static 
head. 

2.  To  determine  if  the  Enger  and  Levy  Equation  can  be  applied  to 
eight-inch  irrigation  pipes. 

3.  To  determine  how  any  empirical  relations  obtained  compare  with 
other  independently  observed  data. 

4.  To  determine  a  design  procedure  for  ensuring  uniform  flow. 
Theoretical  analysis  made  as  part  of  the  investigation  led  to  the 

establishment  of  the  relation 

C,  =  sine/  /~ 2~~ 
d 

The  angle  of  discharge  was  believed  to  be  directly  influenced  by  both 
the  velocity  of  water  approaching  the  orifice  and  the  static  pressure 
at  the  orifice.  However,  no  theoretical  relationship  could  be  established. 
An  experiment  in  pursuit  of  the  above  mentioned  objectives  was  conducted. 
For  the  purpose  of  the  experiment,  an  irrigation  pipe,  30-feet  long 
with  six  circular  holes  drilled  in  it  was  used  as  a  test  pipe.  It  was 
connected  to  another  pipe  as  part  of  a  simulated  pipe  irrigation  system. 
Water  was  run  through  the  system  with  the  aid  of  a  centrifugal  pump.  The 
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discharge  out  of  each  orifice,  within  a  certain  recorded  time,  was  col- 
lected in  a  can  and  weighed.  The  head,  again  at  each  outlet,  was  read 
from  a  set  of  manometers  connected  to  pressure  tappings  at  the  various 
orifices.  Four  runs  were  made  with  each  of  the  two  orifice  sizes  used  in 
the  experiment.  The  flow  rate  was  different  for  every   run. 

The  data  collected  in  the  experiment  was  used  to  determine  the  dis- 
charge coefficients  and  the  approach  velocities  towards  each  orifice. 
The  data,  also,  helped  confirm  the  theoretical  suggestion  that  the  dis- 
charge coefficient  of  an  orifice  is  influenced  by  both  the  head  at  and 
the  approach  velocity  towards  it.  It  further  revealed  the  fact  that  the 
relation  between  the  coefficient  and  these  influencing  factors  is  not  a 
linear  one.  Statistical  analysis  showed  that  a  cubic  model  is  the  best 
fit  for  the  data  in  establishing  a  relation  between  the  coefficient  and 
the  two  quantities,  pressure  head  and  approach  velocity.  When  it  was 
carried  further,  the  analysis  precipitated  two  prediction  equations.  One 
gave  the  discharge  coefficient  as  a  function  of  both  the  pressure  head 
and  the  approach  velocity  towards  the  orifice.  The  other  expressed  the 
coefficient  as  a  function  of  only  the  velocity  of  approach.  The  two 
equations  were  expressed  as: 

Cd  ■  0.5836  +  0.3723V  -  0.01098hV  -  0.346V2  +  0.1084V3 


and 


Cd  =  0.5883  +  0.3106V  -  0.3141V2  +  0.0898V3 


Both  were  found  to  predict  equally  well  within  the  range  of  values  of  the 
static  heads  and  the  pipe  velocities  observed  in  the  experiment. 

The  heads  at  each  outlet  were  established  by  series  of  calculations. 
Hence,  the  difference  between  the  required  and  the  calculated     heads 
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were     computed.    The  orifices  were  first  taken  to  be  at  the  same  ele- 
vation with  respect  to  an  arbitrary  datum.  The  differences  in  head  then 
was  translated  into  an  elevation  difference.  The  slope  of  the  line  of 
least  square  through  the  points  corresponding  to  the  elevation  at  the 
various  orifices  was  taken  to  be  the  optimum  slope.  It  was  found  that 
Enger  and  Levy's  Equation  could  not  be  used  with  this  methodology.  A 
computer  program  was  written  for  easy  application  of  it. 
It  was,  therefore,  concluded  that: 

1.  The  coefficient  of  discharge  could  be  expressed  as  a  function 
of  either  the  approach  velocity  only  or  both  the  approach  velocity 
and  the  static  head. 

2.  The  Enger  and  Levy  Equation  can  be  applicable  to  the  case  of 
eight-inch  pipe  if  good  assumptions  can  be  made  for  the  coefficient 
at  the  dead  end. 

3.  The  prediction  equations  obtained  as  a  result  of  this  study 
predicted  satisfactory   values  of  orifice  discharge  for  different 
conditions  of  flow  and  different  pipe  and  orifice  sizes. 

4.  A  methodology  has  been  established  for  determining  the  optimum 
slope  needed  to  provide  a  uniform  flow,  for  all  practical  purposes, 
out  of  various  orifices  of  a  gated  pipe. 
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SUGGESTIONS  FOR  FUTURE  RESEARCH 

This  study  did  not  include  a  wide  range  of  static  pressures  at  the 
different  orifices.  Therefore,  the  full  extent  of  the  interplay  between 
a  static  regain  due  to  momentum  decrease  and  the  possible  loss  of  pressure 
head  due  to  branching  of  the  flow  could  not  be  determined.  Hence,  a 
further  investigation  to  determine  this  shall  help  eliminate  some  of  the 
assumptions  being  made  in  calculating  the  static  head  at  various  points 
on  the  pipe,  thus  making  the  process  of  determining  the  optimum  slope 
more  accurate. 

The  pipe  used  in  this  experiment  had  plane  circular  orifices  cut 
with  a  drill  press  while  the  conventional  gated  irrigation  pipes  have 
adjustable  gates,  of  varying  shape,  with  some  rubber  fittings  on  the  in- 
side. How  much  this  difference  in  the  nature  and  shape  of  the  orifices 
can  affect  the  applicability  of  the  results  of  this  study  calls  for  a 
further  investigation. 

The  ratio  of  orifice  to  pipe  diameters  has  been  suggested  as  a  pos- 
sible contributor  to  the  variation  in  the  discharge  coefficient  of  ori- 
fices. Experimental  study  can  provide  evidence  for  either  the  justifi- 
cation or  the  dismissal  of  any  assumptions  that  might  have  been  made  with 
respect  to  this  ratio.  Thus,  it  may  contribute  to  a  better  understanding 
of  the  problem. 

Theoretical  analysis  made,  as  shown  earlier  herein,  has  indicated 
a  direct  relation  between  the  discharge  coefficient  and  the  angle  of 
exit  of  the  flow  out  of  the  orifice.  A  further  research  to  determine 
the  factors  affecting  it  are  essential. 
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Finally,  the  range  of  the  velocities  observed  in  this  experiment, 
due  to  the  constraints  of  available  facilities,  was  relatively  narrow 
as  compared  to  those  obtained  in  the  field.  So,  a  research  project,  per- 
haps in  the  field,  that  can  include  the  range  of  flows  for  irrigation 
systems  shall  help  determine  whether  or  not  the  prediction  equation 
obtained  in  this  study  has  to  be  confined  to  the  range  of  values  of 
velocity  observed. 
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APPENDIX  A 


Computer  Program  Developed 
for 
The  Determination  of  the  Optimum  Slope 

that  Would  Provide  Uniform  Flow 
Out  of  the  Orifices  of  a  Multi-Outlet 
Irrigation  Pipe. 
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THIS    CCfFUTER    PRrCFAHME    WAS    WRITTEN    5 Y    SAlIHj    S.    A2UEAXAR 
TC    DEMONSTRATE    FTV.    =  CR    J    PAMFCIC     IRRIGATION    PIFEtAN 
OPTIKUf    StCFE    THAT    WOULD    *CST    CICSELY    APPRCMMTE    EC  L/L 
■DISCHARGE    Tt-RCLGt-     IIS    VARJCLS    OF  I  F I C E S  .  COLL C    EE    C£TEFI*iNEC 


DEFINITION    CF    VAR1AFLES    LSE  C     IN    THIS    PRCGRA^E: 

TCTFLC  =  TI-E    FATE    CF    HCW    CF    WATER     IMSICE    A    PIPE    hITMh    A    SEGfEM. 

TGTAI_C  =  TI-E    TOTAL     FLCV.    ENTERING     Tl-E    PIPE. 

INIHEC=TFE    INITIAL     STMIC    FEAC    -LST    SEFCRE    TFE    FIRST    CP  IP  I C  £ . 

TCTFED  =  Ti-E    TOTAL     F  E  A  C    ECTH    CYNAMC     ANC    STATIC     U    ANY    PC1NT. 

Pf-EC8E!  I)«TFE    STATIC    i-EAC    JlST     5EFOFE     AN    CRIFICE     ill.' 

Pt-ECAF  (  I  )  =  T  5-  E    STATIC    HEAD    JLST    AFTER    AN    ORIFICE     (I) 

PVEAH  =  ThE  .  STATIC    l-EAC     AT     AN    CRIFICE. 

CK  =  hAZEN-WULIAHS'     FR1CT1CN    FACTCF. 

DFIF£=TI-E    CIAVSTEF    Cf    Tl-E    PIPE. 

ORFC*T>E    CIAWETFR    CF    THE    ORIFICE. 

A  =  THE    CROSS-SECTIONAL    AREA    CF    Tt-5    PIPE. 

P=ThE    FEPIfETEP    CF    T(-E    PIPE. 

R=THE    I-YCRALLIC    PACKS    CF     Tl-E    PIPE. 

AC=ThE    AREA    CF    Tt-E    CF1FICE. 

DCCEFF^ThE    CISCI-APGE    COEFFICIENT    CF    THE    ORIFICE. 

Q(I)  =  Tt-E    CISCMPC-E    CUT    CF    AN    CRIFICE. 

V«THE    VELOCITY    CF    AFFRCACI-    TCVARCS    AN    CRIFICE. 

VFEAO  =  VELCC  ITY    h  E  A  0    WITHIN    A    SECTION. 

L  =  LENGTI-    CF    SPACING     EETwEEN    CFIflCES. 

FRCLCS  =  TI-E    FRICTICN    LCSS    BETWEEN    TWC    CGNSECLTIVE    CRIFICES. 

AVEHEC  =  AVERAGE    CF    t-EACS    BEFORE     ANC    AFTER    AN    CRIFICE 

PERC  =  PEPCENTAGE     CEVIATICN     FFCf    Tl-E    AVERAGE    FEAC 

SUI»CSC=SL*     CF    SCLARES 

ELD1FF=C1FFEP-ENCE    BETWEEN    RESLIREO    ANC    ACTUAL    l-EAC    AT    AN    CRIFICE 

B=THE    SLCPE    CF    THE    LINE    CF    LEAST     SOLARES 
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Table  6.  Comparison  of  Calculated,  Required  and 
Adjusted  Heads  and  Discharges. 
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APPENDIX  B 


Sample  Calculations  for  Tables  1  and  2 
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Sample  Calculations  for  Tables  1  and  2 

In  Table  1  run  number  1,  the  weight  of  the  sample  discharge  out 
of  the  orifice  no.  1  at  the  dead  end  was  measured  as  201.75  pounds  and 
was  collected  in  179.9  sees.  Therefore: 

"  =    m.lMll   x  wj  ibs/ft3   x   448-83  a»/cfs 

or 

q  =  8.07  gpm 
where: 

q  =  the  discharge  out  of  the  orifice 
The  approach  velocity  towards  the  orifice  no.  1  is: 

»■  \ 

where: 

V  =  approach  velocity  toward  the  orifice 

Q  =  the  total  discharge  inside  the  pipe  approaching  the  orifice 

A  =  the  cross-sectional  area  of  the  pipe 

Therefore: 

V  =  1    =     8.07  gpm/(448.83  gpm/efs) 
A       v   (8/12  ft)2/4 

or 

V  =  0.051479  cfs 
The  discharge  coefficient  is  given  as 


Cd  =  q/(a/2gh  ) 
where: 

C  .  =  coefficient  of  discharge 

a     =  cross-sectional  area  of  the  orifice 

g     =  acceleration  due  to  gravity 

h    =  the  head  at  the  orifice 
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Therefore: 


8.07  gpm/(448.83  gpm/cfs) 


77(13/(16  X  12))74   •   /  2  X  32-2  X  1-2683 


or 


Cd  =  0.5522 
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ABSTRACT 

Adjustments  have  to  be  made  manually  in  order  to  approximate  equal 
flew  out  of  the  various  gates  in  gated  pipe.  Complete  automation  of  the 
pipe  system  can  not  be  attained  until  a  way  can  be  found  to  solve  this 
problem. 

This  study  was  undertaken  to  solve  the  problem  of  non-uniform  flow 
through  gates.  Studies  made  by  many  people  in  a  similar  or  related  area 
and  well  established  theoretical  analyses  have  indicated  the  dependence  of 
the  discharge  out  of  any  orifice  upon  the  static  head  at  the  orifice  and 
its  discharge  coefficient.  Earlier  studies,  also, suggested  the  effect 
of  approach  velocity  toward  an  orifice  and  the  head  at  the  orifice  upon 
the  discharge  coefficient.  There  was  only  one  empirical  formula  given 
in  the  literature  which  related  the  coefficient  with  the  dynamic  and 
the  static  heads.  The  shortcoming  of  this  relationship  lies  in  the  fact 
that  the  coefficient  of  discharge  for  the  orifice  at  the  dead  end  must  be 
assumed  before  the  coefficient  for  other  orifices  can  be  determined. 

Investigation  was,  therefore,  made  with  the  objectives  of  determin- 
ing the  relationship  between  the  orifice  discharge  coefficient  and  the 
approach  velocity  towards  the  orifice  and  between  the  coefficient  and  both 
the  approach  velocity  and  the  static  head  at  the  orifice.  Further,  the 
applicability  of  the  only  empirical  relationship  available  from  the  lit- 
erature on  the  conditions  of  flow  with  an  eight-inch  irrigation  pipe 
was  investigated.  The  accuracy  of  predictions  by  the  empirical  relation- 
ships obtained  as  a  result  of  this  study  were  also  to  be  determined. 
Finally,  a  procedure  for  using  the  relationships  obtained  in  designing 
optimum  slopes  of  pipes  to  provide  uniform  flow  was  sought. 


A  simulation  of  field  conditions  was  carried  out  in  the  laboratory. 
An  eight-inch  pipe,  thirty  feet  long  was  used  as  the  test  pipe.  The 
friction  factor  of  the  pipe  was  determined  by  running  water  through  and 
determining  the  head  loss  between  two  points,  one  near  the  inlet  end  and 
the  other  near  the  outlet  end.  Circular  holes  were  then  drilled  in  the 
pipe  to  represent  the  gates  that  are  found  in  gated  irrigation  pipes. 
The  pipe  was  connected  to  a  water  source  at  one  end  while  the  other  was 
plugged.  Water  was  then  run  through  such  that  it  was  discharged  through 
the  orifices.  Samples  were  caught  and  weighed  within  some  recorded  time 
intervals  at  each  orifice.  Thus,  the  rate  of  discharge  out  of  each  was 
determined.  Manometer  taps  positioned  opposite  each  orifice  were  used 
in  determining  the  static  head  at  each  outlet.  From  the  data  collected, 
the  approach  velocity  and  the  discharge  coefficient  at  ewery   orifice  was 
calculated. 

Statistical  methods  employing  a  "SAS"  (Statistical  Analysis  System) 
computer  program  was  used  in  analyzing  the  data.  Using  several  statis- 
tical procedures,  two  prediction  equations  were  obtained;  one  expressing 
the  orifice  discharge  coefficient  as  a  function  of  approach  velocity 
while  the  other  expressed  the  coefficient  as  a  function  of  both  the 
approach  velocity  and  the  static  head  at  the  orifice. 

The  two  equations  were  used  to  establish  a  methodology  for  designing 
an  optimum  slope  that  would  give  flows  with  minimum  deviation  from  equality 
of  flow. 


